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One of the fundamental principles of statistical physics is that only partial information about 
a system’s state is required for its macroscopic description. This is not only true for thermal 
ensembles, but also for the unconventional ensemble, known as Generalized Gibbs Ensemble (GGE), 
that is expected to describe the relaxation of integrable systems after a quantum quench. By 
analytically studying the quench dynamics in a prototypical one-dimensional critical model, the 
massless free bosonic field theory, we find evidence of a novel type of equilibration characterized by 
the preservation of an enormous amount of memory of the initial state that is accessible by local 
measurements. In particular, we show that the equilibration retains memory of non-Gaussian initial 
correlations, in contrast to the case of massive free evolution which erases all such memory. The 
GGE in its standard form, being a Gaussian ensemble, fails to predict correctly the equilibrium 
values of local observables, unless the initial state is Gaussian itself. Our findings show that the 
equilibration of a broad class of quenches whose evolution is described by Luttinger liquid theory 
with an initial state that is non-Gaussian in terms of the bosonic field, is not correctly captured 
by the corresponding bosonic GGE, raising doubts about the validity of the latter in general one¬ 
dimensional gapless integrable systems such as the Lieb-Liniger model. We also propose that the 
same experiment by which the GGE was recently observed [Langen et al., Science 348 (2015) 207-211] 
can also be used to observe its failure, simply by starting from a non-Gaussian initial state. 


Introduction - Understanding the physics of quantum 
many-body systems out of equilibrium is one of the most 
challenging open problems today [1, 2], Of central in¬ 
terest is the problem of quantum quenches i.e. abrupt 
changes of the Hamiltonian parameters of a closed quan¬ 
tum system [3], especially in one-dimensional (Id) inte¬ 
grable systems where the study of quantum dynamics 
has led to intriguing discoveries, like the experimental 
observation of lack of thermalization [4] and the theoret¬ 
ical prediction of the Generalized Gibbs Ensemble (GGE) 
[5] which has recently been observed experimentally [6]. 
The GGE is expected to describe the equilibration of lo¬ 
cal observables after a quantum quench in an integrable 
system by taking into account all constraints associated 
to its conserved charges [3, 5, 7-45]. While in its stan¬ 
dard form it was constructed exclusively out of local 
conserved charges [28, 29], it has been recently shown 
that quasi-local charges must also be included [46-60]. 
This gives an answer to the fundamental question of how 
much information about the initial state survives in the 
final values of local observables. Other aspects of in¬ 
terest are the asymptotic behavior towards equilibrium 
[13, 30, 34, 41, 61] or the restoration of symmetries of 
the post-quench Hamiltonian that are absent in the ini¬ 
tial state [62, 63]. 


solely in terms of the initial two-point function. This 
proves the validity of the GGE which for free systems is 
a Gaussian ensemble containing precisely the same infor¬ 
mation. The argument hinges upon the cluster decompo¬ 
sition property of the initial state and the presence of a 
non-vanishing post-quench mass. In the present work we 
study the case of massless evolution. We find that, unlike 
the massive case, the large time limit of physical observ¬ 
ables retains memory of non-Gaussian initial correlations 
in the form of spatial averages of multi-point connected 
correlation functions (CCFs) of all orders. Therefore the 
equilibrium cannot be described by a Gaussian ensem¬ 
ble and the GGE fails, unless the initial state is Gaus¬ 
sian itself. Additionally, we derive exact results for the 
asymptotic behavior of physical observables and correla¬ 
tion functions at large times and distances. We find that 
symmetries that characterize the post-quench Hamilto¬ 
nian but are in general absent in the initial state are 
restored by the evolution. Lastly we compare our con¬ 
clusions with earlier studies discussing the scope of their 
validity and proposing an experimental implementation. 

Model and quench protocol - We consider a quantum 
quench in a Id bosonic system, starting from any massive 
field theory, free or interacting, to the free massless field 
theory. Thus the post-quench Hamiltonian is 


In [34] it was shown that in the case of a quantum 
quench from an interacting to a free massive bosonic 
field theory in Id, the evolution eliminates memory about 
non-Gaussian correlations in the initial state so that the 
large time values of any local observable can be expressed 


H = 


dx 




(1) 


while the pre-quench one is Hq = H + f dx V{4>(x)), 
with some local ‘potential energy’. The initial state 
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is the ground state of Hq and we fix the field <j> so that its 
initial expectation value vanishes (0) o = 0 (the subscript 
‘0’ denotes the initial state). We also assume that the 
Hamiltonian H 0 is massive, meaning that its ground state 
exhibits short-range correlations i.e. all connected corre¬ 
lation functions (CCFs) decay exponentially with the dis¬ 
tance within a finite correlation length fo = 1 /too where 
mo is the (renormalized) mass of its lightest particle [64] . 

The free massless Hamiltonian (1) describes the sim¬ 
plest CFT, that is, the simplest Id critical model [65]. Its 
ground state exhibits logarithmic 0 correlation functions 
(0(0,O)0(x, t)) gs = -log(|a;| 2 -|t| 2 )/(47r)-l-const. Physi¬ 
cal observables are however expressed in terms of the so- 
called ‘vertex operators’ V a {x,t ) = e la 'K x ’ t ) or derivatives 
of </>, which in the ground state of H exhibit power-law 
decaying correlation functions, as expected for a critical 
model, e.g. (V a (0,0)V- a (x,t)) gs ~ (|x| 2 - |t| 2 )-“ 2 /( 47r ). 
All other ground state multi-point functions can be de¬ 
rived from the above using Wick’s theorem, since the 
ground state of any free Hamiltonian like H is Gaussian. 
Wick’s theorem for vertex operators reads 


_ g-| I2i a i <^i> 


( 2 ) 


Excitations of H correspond to massless particles with 
linear dispersion relation £/ = |fc| (to simplify notation, 
we have set the speed of sound v = 1). 

Note that H is symmetric under the continuous trans¬ 
formation 0 —> 0 + e. This is reflected in the fact that 
correlation functions of vertex operators vanish unless 
they satisfy the ‘neutrality condition’, more explicitly 
(TL v ai(xi,ti)) gs = 0 for J2< a i 7 ^ In particular one- 
point functions with a ^ 0 vanish. The initial state in¬ 
stead, being the ground state of an arbitrary Hamilto¬ 
nian, breaks in general this symmetry. Similarly, H is 
symmetric under the discrete transformation 0 —>- — q i. 

Evolution of observables - We are interested in the 
time evolution of physical observables, in particular equal 
time correlation functions of vertex operators at different 
points, focusing on their asymptotic form at large times 
t after the quench. The evolution of the one-point func¬ 
tion can be derived most conveniently from the cumulant 
expansion 


(V a (x,t)) = exp 


DHVm,, 

; n\ 


( 3 ) 


where the double brackets ((...)) denote connected cor¬ 
relation functions. We work in the Heisenberg picture, 
i.e. all correlation functions (...) refer to initial state 
expectation values of time-evolved operators. Similarly 
the two-point correlation function of vertex operators at 


different points x\ ^ X 2 is given by 


(V a (x 1 ,t)V b (x 2 ,t)) 


= exp 


w - n 

(W’fci a) +h(x2, *))"} 


( 4 ) 


The relation between CCFs and standard ones is that 
between joint cumulants and moments. 

From the above we see that it is sufficient to study 
the evolution of each CCF of 0. The equation of motion 
for 0 after the quench is simply the Id wave equation in 
infinite space 

4>(x,t) = d 2 x cj){x,t) , t> 0 (5) 

with two initial conditions 0(x, 0) and 0(x,O). The solu¬ 
tion is given by d’Alembert’s formula 

1 1 r x+t 

0(x, t) = — (0(x — f, 0) + 0(x +1, 0)) + — / dx' 4>(x', 0) 

2 2 J x _ t 

( 6 ) 

or in terms of the Green’s function 

0(x,f) = ^ f dx' G ( - S \x — x',t)f)^ 1 ~ s \x',0) (7) 

s-0,1 2 


where G(x, t) is the (retarded) Green’s function for the 
Id wave equation in infinite space, that is G(x , t) = 

f dk e lkx (sin kt)/k = 40(t — |x|) (we use the nota¬ 
tion / (0) (i) = f(t),f W {t) = f(t) = d t f(t)). Note the 
light-cone structure of G(x,t)- It then follows that the 
time evolution of CCFs is 

j / Tt \ \ a Tl TL 

n#^) = / n yi 

\\^_l // J 2=1 i— 1 s i —0 1 

x (n^- o w.o)) ( 8 ) 

It is now clear that using the cumulant expansion and 
the last formula, we can directly express any correlation 
function of vertex operators in terms of a suitable convo¬ 
lution of initial CCFs of 0 and the large time asymptotics 
is controlled by the scaling of the above convolution, that 
is, by the large distance behavior of initial correlations. 
Below we study separately the case of Gaussian and non- 
Gaussian initial states. 

Gaussian initial states - In this special case, all n- 
point CCFs of 0 with n > 3 vanish. This holds when 
H 0 is essentially non-interacting (i.e. either explicitly 
free, as when V(0) = 4 toq 0 2 , or interacting that can 
be mapped into free). Since we have set (0(x, 0)) = 0 we 
now have ((0(0,0)0(x, 0))) = (0(0, O)0(x, 0)), i.e. the only 
non-vanishing initial correlation function is the two-point 
function. 

The one-point function of vertex operators is 
(14(0, t)) = e-5 Q2 ^ 2 (°’‘)> and from (6) and (8), tak¬ 
ing into account the exponential large-distance de¬ 
cay of initial correlation functions, we find that the 
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leading large-time behavior is given by (4> 2 (0,t)) = 
f dr(<p( 0, 0 )4>(r, 0)) + ... where the dots denote terms 
that increase slower with time for t £o- The last inte¬ 
gral tends to a finite value A2 as t —> 00, since the initial 
two-point function (</>(0, 0)(f>{r, 0)) decays exponentially 
with the distance r. Therefore the leading large time 
behavior of the vertex operator one-point function is 

(14(0,f)} oce-^° 2A2t (9) 

where, using time translation invariance for t < 0 and 
the pre-quench equation of motion <p(x,t) = d 2 <p(x,t) — 
V'(<f>(x,t)) evaluated at t —> 0 _ , we obtain 

/ +00 _ _ r+co 

dr(<j>{ 0,0)^(r,0)>= / dr(^(0, 0)V'(^(r, 0))) 

-OO J — OO 

( 10 ) 

We therefore find that the one-point function decays 
exponentially with a decay rate A2 determined for any 
Gaussian initial state by (10). Note that in the contin¬ 
uum limit A2 is typically an ultraviolet-divergent quantity 
[3] meaning that it depends on microscopic details of the 
actual physical system (lattice spacing etc.). 

Similarly we calculate the two-point function 
of vertex operators, which is ^ e la 0( o,t) e ib</>(x,t)^ _ 
exp [—|(a + b) 2 (cf> 2 (0, t)) - ab{{<j>(0,t)<j>(x,t)) - (0 2 (O, i)>)] 
For large times t^> £o,t > |x|/2 we find t)<j>[x, t)) — 
(0 2 (O,i)) = -|/d a| dxi /(j x| dx 2 (<j)( 0 , 0 )(j)(x 2 -x 1 , 0 ))+..., 
which for |x| £0 tends to — ^A2 |a?|- Essentially the 

asymptotics for \x\,t 3> £0 can be easily calculated 
substituting 0)<^(x2,0)} —> A2<5(x2 — x\). From 

the above we finally find 

(V a (0 : t)V b {x,t)) oce-^+^t+iabA^zl (11) 

i.e. two-point functions of vertex operators decay expo¬ 
nentially with time, unless b = —a in which case they 
equilibrate to nonzero values. The large time two-point 
function (14(0)E_ a (a;)) t _ >oo oc e _ « a A2 I X I decays expo¬ 
nentially with the distance with a decay rate controlled 
by the same parameter A2 as the time decay. 

For multi-point correlation functions we similarly find 

( n *)) « e"3 A2( ^ a d 2( +i A 2 Z i<3 i-Xil 

( 12 ) 

in agreement with Wick’s theorem (2) which holds be¬ 
cause both the initial state and its evolution are Gaus¬ 
sian. The above shows that all multi-point correlation 
functions of vertex operators decay exponentially with 
time, unless they satisfy the neutrality condition, that 
is, the massless evolution restores the symmetry under 
(f> —> <f> + e transformations. 

Non-Gaussian initial states - In the more general 
case of a non-Gaussian initial state, like the ground 
state of a genuinely interacting Hamiltonian Ho (non¬ 
parabolic potential V{<j>) = \m g</> 2 + c n<t> n /n !) 


the large time asymptotics of vertex correlation func¬ 
tions can be calculated following the same arguments 
as above, taking into account that for a massive self¬ 
interacting pre-quench Hamiltonian, all initial CCFs de¬ 
cay exponentially within a range of the order of the cor¬ 
relation length £ 0 . The vertex operator one-point func¬ 
tion is given by (3) with <f> correlation functions given 
by (8) where at large times we now find (0™(O,t)J> = 

2 1 -”inr=2/- 4 r^^(0,°)nr=2^ ) 0))) + ... and, as 
before, for f > ( 0 we can replace (n"=i 0))) 

A — Xi) where 

r + °° n //. n w 

x n = nM*< o ’°>nto>°>> ( i3 ) 

i =2 i=2 


Thus we finally find 


(K(0 ,t)) oc exp 



(-»") 

u!2 n_1 



t 


(14) 


The expression in round brackets is a series of spa¬ 
tial averages of all initial multi-point CCFs. Since 
all A„ are real numbers, it is generally a complex 
number, whose real/imaginary part is the sum of 
even/odd terms respectively. This means that the 
one-point function decays exponentially with time with 
decay rate 7(a) = E)Ei( _1 )" +la2 "A2n/[(2n)!2 2Ti_1 ] 
and also exhibits oscillations with frequency to (a) = 
Er=i(- 1 ) na2 ” +1 A2n+i/[(2n + l)!2 2 "]. If the pre- 
quench Hamiltonian is symmetric under reflections (j) —>- 
— (j>, then all odd correlation functions vanish and the 
evolution is purely decaying. The quantities X n clearly 
depend on all parameters of H$\ e.g. if perturbation the¬ 
ory describes correctly its ground state, then at first order 
in the couplings c n they are A2 cx mo and A„>3 oc c n . 

Similarly for the two-point function, from (4) and (8) 
we find 


(V a (0,t)V b (x,t)) oc exp - ( — y^(o + b) r 


-Xn t 


-(EKa+sr-K+ni^MM 


(15) 


which decays with time, unless b = —a in which case it 
equilibrates to (I4(0)H_ a (a;)) t ^ oo oc e _7 ^ a ^ a L Note that 
the equilibrium values depend not only on the initial two- 
point CCF of (f) as for Gaussian initial states, but also on 
all even higher ones, through the values of the indepen¬ 
dent parameters X 2 n - This shows that information about 
non-Gaussian features of initial correlations survive at 
t —» 00. It should also be noted that the evolution re¬ 
stores the symmetries under both <j) —> <fi+e and </> —> —(j). 

Failure of the bosonic Gaussian GGE - The GGE den¬ 
sity matrix in one of its standard forms for free systems, 
is constructed using as integrals of motion the number 
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operators n k , i.e. it is written as p g ge oc e~^ k ^ knk with 
fixed by the condition that the values of the charges 
in the GGE equal their initial values {n/ c ) gge = (n k )o [3] . 
The operators n k are linear combinations of the local in¬ 
tegrals of motion and, in any finite lattice system, equal 
in number to them [28, 29]. Written in this form, de¬ 
spite being non-local, the GGE takes into account the 
quasilocal charges of continuous models defined in [59] . 

Since this is a Gaussian density matrix, the GGE 
prediction for any multi-point correlation function of 
vertex operators satisfying the neutrality condition 
is given by Wick’s theorem (2) (fl,e ,ail ^ x< )) gge = 

with ]T\ai = 0 (as 
already mentioned, if the neutrality condition is not 
satisfied the correlation functions vanish at t —> oo). 
This means that all information about the GGE is 
included in its two-point function of <f >, which can 
be expressed in terms of (n k )o as (</>(0)</>(x))g ge = 

^ J dke lkx [l +((n k + n- k )) 0 \/{2E k ) where E k = |fc|. 

In terms of initial correlation functions, using the pre¬ 
quench equation of motion, we find 

(</>(O)0(x)) gg e = {(j>(0)^(x)) o - J dr ^ 4 7 ’^ (O)V , (0(r))) o 

(16) 

which can be seen as a GGE version of the virial theo¬ 
rem that relates the temperature of a thermal ensemble 
with the inter-particle potential energy. The asymptotic 
form at large |x| is (^>(0)cj)(x)} gg e ~ —\^ 2 \x\ + const., 
where A2 is given by (10). Note that in the special case 
V(</>) = \m^(j) 2 which corresponds to a free massive pre¬ 
quench Hamiltonian, the parameter A2 equals mo/2 and 
the GGE is in agreement with our t —> 00 results for 
Gaussian initial states. 

As already mentioned, the above ensemble is by con¬ 
struction Gaussian and captures only information about 
the initial two-point function of (j> but not of higher or¬ 
der initial CCFs. Therefore its predictions do not give 
the correct t —> 00 values of vertex operators. This is 
also true for the generalization of the GGE that includes 
the quasi-local charges of [59] since that too is Gaussian. 
The fact that local observables retain memory of all even 
initial CCFs suggests that an ensemble that describes cor¬ 
rectly the t —> 00 limit would also include all products of 
n k as conserved quantities 

( OO t \ 

Ai,..,^II n b ( 17 ) 

j —1 / 

This is clearly equivalent to the Diagonal Ensemble, i.e. 
the expansion in terms of projections onto all energy 
eigenstates [7]. This ensemble includes all information 
about the initial state that contributes to the infinite 
time average of any observable in a general system [2, 7] 
and is therefore both non-economic and impractical. It is 


still possible that a suitable non-Gaussian generalization 
of p gge that would include extra conserved charges that 
are special of CFT can describe correctly the large time 
limit. In either case, the free massless Id field theory 
apparently retains the maximum possible memory of its 
initial state. 

Note that, while it has been very early pointed out 
that the GGE does not capture information about initial 
correlations between different n k [16, 20, 66], in the ther¬ 
modynamic and large time limit local observables typi¬ 
cally lose memory of such correlations [16, 67], which is 
the deeper reason that allows their economic description 
through the GGE. In contrast, in the present problem, 
information about all such correlations survives in the 
above limit and remains accessible through local mea¬ 
surements. To the best of our knowledge this is the 
first example of its kind. An intuitive explanation of 
this memory effect is that, due to the ballistic nature of 
the evolution, entangled clusters of quasiparticles carry 
the initial correlations intact up to infinite time. 

Discussion of results - The massless Gaussian field the¬ 
ory described by (1) is a prototypical model of Id phys¬ 
ical systems with gapless (phononic) excitations. Indeed 
the ‘bosonization’ method shows that the low-energy be¬ 
havior of many interacting systems, both fermionic or 
bosonic, can be described by the harmonic-fluid or Lut- 
tinger liquid (LL) theory [68-70]. These include the 
Lieb-Liniger model (which is supposed to describe the 
Quantum Newton’s Cradle experiment [4]) and the gap¬ 
less phase of many lattice models (e.g. Bose-Hubbard 
model in the superfluid phase) and spin chains (e.g. 
the gapless phase of the XXZ model) in the scaling 
limit. The Hamiltonian of the Luttinger model is Hll = 
\v f dx (^-(d x 0) 2 + -^n 2 ) where the density n(x) and 
the phase field 9{x) obey canonical commutation rela¬ 
tions and the parameters v and K are the sound velocity 
and Luttinger parameter respectively. The above form is 
equivalent to (1) as can be seen by rescaling the field and 
the space-time coordinates. 

In this context our findings show that the GGE as in¬ 
terpreted in the LL approximation, i.e. constructed out 
of the local charges of the bosonic density-phase fields, 
does not describe correctly the steady state correspond¬ 
ing to a general initial state. We emphasize that LL 
theory is only a low-energy approximation of the original 
models, while a quantum quench may create arbitrarily 
high-energy excitations. Therefore in this context our re¬ 
sults do not necessarily mean that the GGE conjecture 
itself is incorrect: it is possible that the GGE constructed 
out of the exact conserved charges of the original model 
is correct. In particular in the case of an instantaneous 
interaction quench in the Lieb-Liniger model, the pres¬ 
ence of high-energy excitations [31, 33] makes the LL ap¬ 
proximation non-applicable and at least the special cases 
studied in the literature seem to suggest that the GGE 
expressed in terms of the exact Lieb-Liniger conserved 
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charges is correct [33, 34, 38, 44, 45, 71]. Neverthe¬ 
less the LL approach turns out to describe correctly the 
quench dynamics in several other cases [72-74] . Quantum 
quenches in LL have been studied in [15-17, 75, 76] and 
a detailed comparison of our results with earlier works is 
presented in the Supplemental Material. 

Note that our method [34, 61] for the derivation of 
quench dynamics has certain advantages in comparison 
with others. By exploiting the connection between ini¬ 
tial and evolved fields, we link the large time asymptotics 
of correlations directly to the long distance asymptotics 
of initial correlations, which are correctly described by 
Renormalization Group theory. Even though the appli¬ 
cability of this approach seems to be restricted to models 
with Gaussian evolution, in fact it is based on fundamen¬ 
tal properties of quantum field theory that are also valid 
in the presence of interactions, more precisely the exis¬ 
tence of a suitable set of local fields whose time evolution 
can be expressed as a convolution of initial local fields. 

In particular the conclusion that the GGE fails in the 
case of a genuinely interacting to free massless quench in 
Id relies on the ballistic form of the evolution (6) which 
is linked to the gaplessness of the energy spectrum: if 
the post-quench Hamiltonian was gapped instead, the 
Green’s function in (7) would decay with the distance in¬ 
side the light-cone and the quench dynamics would lead 
to the GGE [34]. Therefore our result is expected to be 
robust under the insertion of irrelevant interactions to the 
post-quench Hamiltonian, as long as they do not affect 
the ballistic form of the evolution. Such perturbations of 
the Luttinger model preserve the gaplessness of the spec¬ 
trum, even though they modify the dispersion relation 
from linear to nonlinear [77]. The same is also expected 
to be true for quenches in Id CFT in which there always 
exist local operators whose evolution is governed by the 
wave equation (5) (e.g. the energy-momentum tensor). 

Experimental implementation - Ironically, our find¬ 
ings about the failure of the GGE can be experimentally 
checked in the same way that the GGE was actually ob¬ 
served in [6], simply by starting from a non-Gaussian 
initial state, i.e. the ground state of a genuinely interact¬ 
ing Hamiltonian. We will show that this is feasible using 
present technology in the same experimental setup [78]. 

One-dimensional bosonic models can be realized in cold 
atom systems [79] and quantum quenches can be imple¬ 
mented by tuning the system’s parameters [80]. The 
most suitable experimental setup for this purpose is a 
system of two coupled Id gases of interacting bosons pre¬ 
pared by condensate splitting [81-85]. The low-energy 
physics of this system can be derived through bosoniza- 
tion [69, 81-84, 86-89]. The density n and phase field 
9 of the antisymmetric modes are described by the sine- 
Gordon model, with the interaction a controlled by the 
tunneling coupling J between the two Bose gases. For 
J = 0 the interaction vanishes and the sine-Gordon 
model reduces to the Luttinger model i.e. the system 


corresponds to free gapless phonons. In the opposite 
limit of large J, solitonic excitations between neighbor¬ 
ing vacua are suppressed and the system is described 
by essentially free gapped excitations above one of the 
degenerate vacua. The rapid tunability of the param¬ 
eters of the experimental system allows the implemen¬ 
tation and efficient study of quantum quenches. The 
subsequent evolution of correlation functions is analyzed 
through time-resolved measurements of matter-wave in¬ 
terference patterns after time-of-flight expansion of the 
gas, averaged over many repetitions [90]. Such measure¬ 
ments provide direct data for the multi-point correlations 
of the local phase difference 9[x) between the two con¬ 
densates. 

These techniques made recently possible the celebrated 
first experimental observation of a GGE [6]. More pre¬ 
cisely it was demonstrated that more than one ‘effective 
temperatures’ are required for the description of the long 
time steady state of the system, which is a clear sign of 
absence of thermalization and equilibration to some type 
of GGE. In the regime of parameters used in the exper¬ 
iment, the system is well-described both before and af¬ 
ter the quench by the Bogoliubov-de Gennes equations, 
which provide a mean-field (that is Gaussian) approxi¬ 
mation [80, 86, 88, 89]. In particular the post-quench 
Hamiltonian corresponds to the gapless Luttinger model. 
According to our findings, equilibration in this case is cor¬ 
rectly described by the standard GGE since the initial 
state is also Gaussian. 

Deviations from the latter would however be mani¬ 
fested if the initial state were non-Gaussian. In fact 
measurements of CCFs precise enough to reveal non- 
Gaussian features are also feasible in the same exper¬ 
imental setup and have been performed in the ground 
state of the system for a wide range of parameter val¬ 
ues [78]. This suggests that an experimental realization 
of a quantum quench of J from an intermediate value 
to zero would make the observation of such deviations 
possible. Indeed, provided that the effective low-energy 
field theory gives a faithful description of the experimen¬ 
tal system and quench protocol, in this case the initial 
state would be the ground state of the genuinely interact¬ 
ing sine-Gordon model and would exhibit non-Gaussian 
correlations. The parameters X n of (13) can then be 
calculated using the methods of Integrable Field Theory 
[91, 92]. Numerical checks can be performed by means 
of the so-called Truncated Conformal Space Approach 
[93, 94] after a careful analysis of finite size effects. De¬ 
tails will be given in a forthcoming publication. Note 
however that eqs.(13), (14) and (15) are more suitable 
for an experimental check: while it has been difficult to 
compare experimental measurements to theoretical pre¬ 
dictions for the sine-Gordon model, it is quite easy to di¬ 
rectly compare measurements on the initial ground state 
and at large-times after the quench. 

Conclusions - We studied a quantum quench from an 
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interacting to a free massless bosonic field theory in Id, 
deriving analytical results for the large time correlations 
and demonstrating that, in contrast to the massive case, 
the standard Gaussian GGE fails: the system retains 
memory of all initial correlations, which is accessible by 
local measurements. 
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SUPPLEMENTAL MATERIAL 


Comparison with earlier work 

Quantum quenches in the context of Luttinger liquids 
have been studied in [15-17, 72-76] and as a special case 
of CFT in [3, 95]. Even though, as already mentioned, 
LL theory is only a low-energy description of Id gapless 
systems, there is evidence that it captures correctly the 
quench dynamics in several cases [72-74, 96, 97] and cer¬ 
tainly whenever high-energy excitations are suppressed, 
like when the quench is not performed instantaneously 
but within a short time interval. All results of [3, 15- 
17] are consistent with relaxation to the GGE, while the 
more general initial states considered in [95] lead to equi¬ 
libration described by a GGE in which non-local charges 
are also expected to be included. The initial states con¬ 
sidered in [3, 15-17], when specialized to our settings, be¬ 
long to the class of Gaussian initial states, therefore these 
results are in agreement with ours. Indeed the Dirichlet 
initial state used in [3] exhibits vanishingly small correla¬ 
tion length £o ~ m q 1 —> 0 in which limit it is Gaussian. 
In this limit, from (10) we obtain A2 ~ mo/2 recovering 


the results of [3|. 

On the other hand, all quantum quenches in the Lut- 
tinger model studied in [15, 16] correspond, after ap¬ 
plying the bosonization mapping, to quenches within a 
Gaussian model. The same is true for the quench studied 
in [17] since the Luther-Emery point of the sine-Gordon 
model is essentially free. The initial state considered 
in [76] instead is a special type of non-Gaussian state 
and it would be interesting to check if relaxation retains 
memory of its non-Gaussianity, which is not however ad¬ 
dressed in that work as the authors compare only with 
the thermal ensemble predictions. Lastly, the general 
quenches studied in [95] using CFT include also non- 
Gaussian perturbations in the initial state. The obser¬ 
vation that memory of each of them survives at infinite 
times and must be included in the GGE in order for it to 
be correct, is consistent with our result for failure of the 
conventional GGE. 

In the literature an equivalent form of the Diagonal En¬ 
semble (17) has appeared in a different context, the Lieb- 
Liniger model, and has been termed Generalized GGE 
[53]. However, while equilibration of a general initial 
state is always described by the Diagonal Ensemble [7], 
quench initial states, being ground states of local Hamil¬ 
tonians, are special and so far there is no other known 
case of a quantum quench for which this ensemble does 
not reduce to some version of the much more economic 
GGE [60], 


In lattice models, it has been proven in [24, 25] that 
a general initial state satisfying mild clustering require¬ 
ments (as quench initial states do) and evolving under a 
gapless Gaussian Hamiltonian relaxes locally to a Gaus¬ 
sian reduced density matrix. However in that context lo¬ 
cal relaxation refers to single site observables [24], while 
for blocks of neighbouring sites [25] the deviations from 
Gaussianity increase with the block size, meaning that 
in the scaling limit where the lattice spacing is taken to 
zero, correlations at finite distances are non-Gaussian. 
Consequently there is no inconsistency with our results. 
Extending those results to free fermion lattice systems, 
it was more recently shown [98] that, apart from cluster¬ 
ing initial correlations, a second condition for Gaussian 
relaxation is that the evolution exhibits sufficiently delo¬ 
calizing transport. This condition is shown to be typical 
in lattice systems, but is obviously not satisfied in the 
continuous bosonic system studied here, which exhibits 
instead ballistic transport. 

An interesting question is whether the so-called 
Quench Action method [13, 35, 46, 48, 99] would give 
correct predictions in the cases considered here. This 
method approaches the problem of quench dynamics in 
a way different than the GGE and aims at describing the 
thermodynamic properties of the initial state of Bethe 
Ansatz solvable systems in an efficient and faithful way. 
For this reason and although there are no analytical cal¬ 
culations of steady state multi-point correlation functions 
so far, it is not expected that this approach would fail. 


